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We present the results of a search for Higgs bosons predicted in two-Higgs-doublet models, in 

the case where the Higgs bosons decay to tau lepton pairs, using 1.8 fb~^ of integrated luminosity 
of pp collisions recorded by the CDF II experiment at the Fermilab Tevatron. Studying the mass 
distribution in events where one or both tau leptons decay leptonically, no evidence for a Higgs 
boson signal is observed. The result is used to infer exclusion hmits in the two-dimensional space 
of tan (3 versus tua (the ratio of the vaccum expectation values of the two Higgs doublets and the 
mass of the pseudoscalar boson, respectively). 
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Understanding the origin of electroweak symmetry 
breaking is one of the central goals of particle physics. 
The Higgs mechanism [l| in the standard model (SM) 
provides a possible explanation, but the calculated mass 
of the Higgs boson suffers from large radiative correc- 
tions. Remedies for this problem such as supersymme- 
try require at least two Higgs doublets Q and result 
in a more complicated Higgs boson sector than that of 
the SM. The minimal supersymmetric standard model 
(MSSM) Q predicts the existence of three neutral Higgs 
bosons. The MSSM is an example of a Type H two- 
Higgs-doublet model (Type H 2HDM) in which there is 
a light scalar h, a heavy scalar H, and a pseudoscalar A. 
The masses of these states are governed mainly by two 
parameters in the theory, usually taken to be tan/3, the 
ratio of the vacuum expectation values of the two Higgs 
doublets, and itia, the mass of the pseudoscalar. 

In pp collisions at 1.96 TeV center of mass energy 
at the Fermilab Tevatron, Type H 2HDM Higgs bosons 
would be predominantly produced by gluon-gluon fusion 
through a b quark loop 5] or by 66 fusion [6|. The cou- 
plings and masses of the Higgs bosons are such that if 
tan /3 is greater than about 20 and tua is smaller (greater) 
than about 125 GeV/c^, one finds that the h (H) and 
A are degenerate in mass to within a few GeV/c^, and 
are produced with a cross section proportional to tan^ f3, 
while the production of H (h) is suppressed. 

These tan^ /3-enhanced production cross sections can 
be in the range 0.1 - 10 pb depending on the Higgs bo- 
son masses, and are orders of magnitude greater than the 
corresponding ones for a SM Higgs and also the more fa- 
miliar associated production modes of a SM Higgs boson 
with a vector boson. The Higgs bosons decay to fermion 
pairs with a partial width proportional to the fermion 
mass squared; thus the decays (/> — > 66 and </> — > t+t" 
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(with (f> — h, A, H) predominate, with the branching ra- 
tio to 66 approximately 90% and the branching ratio to 
T+r- about 9% for tua > 100 GeV/c^. 

This Letter presents the results of a search for the 
production of Higgs bosons in Type II 2HDM such as 
the MSSM, using data collected with the CDF II detec- 
tor at the Fermilab Tevatron pp collider corresponding 
to 1.8 fb~^ of integrated luminosity. Full details of the 
analysis are available elsewhere [7|; this result supersedes 
our previously published result Q, and is similar to the 
search performed by the DO Collaboration [§| . The anal- 
ysis is sensitive to a region of MSSM parameter space 
which is com plem entary to that explored by the LEP 2 
experiments. [10| 

The analysis presented here uses the tau pair decay 
modes, since it is possible to efficiently trigger on and 
reconstruct the leptons in decays of the tau lepton to 
evD or /ii/P. Indeed, despite the 10 x larger branching 
ratio to 66, the search in the tau mode is more sensitive 
because the SM background is much smaller. 

CDF II [ll| is a general-purpose detector with an over- 
all cylindrical geometry surrounding the pp interaction 
region. The three-dimensional trajectories of charged 
particles produced in pp collisions are measured at small 
radii (< 30 cm) using multiple layers of silicon microstrip 
detectors, and at outer radii (> 30 cm) with a multi-wire 
drift chamber. The tracking system is inside a solenoidal 
magnet with uniform 1.4 T magnetic field oriented along 
the beam direction. Outside the solenoid are the elec- 
tromagnetic and hadronic calorimeters, which are seg- 
mented in pseudorapidity (77) and azimuth in a projective 
"tower" geometry [12]. A set of strip and wire chambers 
located at a depth of six radiation lengths aids in iden- 
tifying photons and electrons from the electromagnetic 
shower shape. Muons are identified by a system of drift 
chambers and scintillators placed outside the calorimeter 
steel, which acts as an absorber for hadrons. The inte- 
grated luminosity of the pp collisions is measured using 
Cerenkov luminosity counters [Tsj . 

We seek events with tau pairs where one or both taus 
decay leptonically (excluding e^e^ and /i^/i^ which suf- 
fer from excessive background from Z/7* production). 
These final states are denoted e + t, fj, + t, and e + fi 
(where "t" here means the reconstructed hadronic part 
of a tau decay). Events with a high-pT (8 GeV/c or 
more) e or candidate plus a high-pT charged track (5 
GeV/c or more) or a second e or ^ (4 GeV/c or more) 
are identified using high-speed trigger electronics and are 
recorded for later analysis. The performance of the trig- 
ger and lepton identification algorithms is described in 
detail elsewhere [13, [3 ■ 

The reconstruction of hadronic decays of tau lep- 
tons ^1^ relies on defining tau signal and isolation region 
cones centered around seed tracks having pT > 6 GeV/c; 
we demand one or three charged tracks in the tau cone, 
and include in addition any tt^ candidates. The main way 
to discriminate between hadronic tau lepton decays and 
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TABLE I: Mean expected SM backgrounds and observed 
numbers of selected events in the final sample. The uncer- 
tainties include all systematic efi'ects, some of which are cor- 
related. 



Background 


e + ^ 


e + T 


fl + T 




605 ± 51 


1378 ± 117 


1353 ± 116 


Z/j* e+e" 


1.5 ± 1.2 


70 ± 10 


negi. 




17.9 ± 4.5 


negl. 


107 ± 13 


dibosons 


11.4 ± 3.7 


4.2 ± 2.1 


3.3 ± 1.8 


ti 


9.1 ± 3.3 


4.0 ± 2.1 


3.3 ± 1.9 


W-|-jet, multijet 


57.1 ± 13.5 


467 ± 73 


285 ± 46 


Total 


702 ± 55 


1922 ± 141 


1752 ± 129 


Observed 


726 


1979 


1666 



to demand no additional charged tracks or 7r° candidates 
in an isolation annulus outside the tau signal cone but 
within 30° of the tau seed track. The half-angle of the 
tau signal cone decreases with increasing visible tau en- 
ergy due to the Lorentz boost of the tau lepton, further 
aiding the discrimination of taus from jets. Additional 
suppression of hadronic jets comes from imposing a mass 
requirement on the tau candidate decay products. Elec- 
trons and muons are removed using information from the 
calorimeters and muon detectors. 

To select the e + t and fi + t events we require an 
isolated e or /i with px > 10 GeV/c, and a r with visi- 
ble hadronic decay products with total pr > 15 GeV/c 
(20 GeV/c for three-charged-pion decays). For the e -t- /i 
channel we require one lepton to have px > 10 GeV/c 
and the other to have pr > 6 GeV/c. 

The main SM contributions to the selected event sam- 
ple include Z/'j* t'^t~, and M^-|-jet events where 
W (with £ — e, fj.) and the hadronic jet is misiden- 

tified as a hadronically decaying r. The VF-|-jet events 
are largely removed by requiring that the missing trans- 
verse energy fix not point along the direction opposite 
the momentum of the £+t system. The remaining iy-|-jet 
background, and all other background stemming from 
jets misreconstructed as taus ("fakes") is estimated from 
events recorded with a jet trigger. There are small con- 
tributions from Z — > e"*"e~, Z p,~ , diboson, and tt 
production. 

The acceptances for signal and the non-fake back- 
grounds are estimated from samples of simulated events 
produced by the PYTHIA event generator [l^ with 
CTEQSl fl§\ parton distribution functions. The Higgs bo- 
son widths and masses are those for an MSSM model with 
tan P = 50. Tau decays are simulated by the tauola 
package [l^. A GEANT-based [l^l model simulates the 
interactions of all final state particles in the detector. 

Table U shows the mean expected contributions of SM 
sources, and the number of observed events in the three 
channels. The uncertainties listed include all systematic 
effects discussed below, including correlations. 

To discriminate a Higgs boson signal from the back- 
grounds, we perform a binned likelihood fit of background 
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FIG. 1: Observed and predicted distributions of rrivis for the 
e + T and fi + t channels (a) and e + /i channel (b). The 
predicted signal distribution (for (j> — h/A/H) corresponds to 
that for niA ~ 140 GeV/c^ signal assuming a value of the 
cross section excluded at 95% C.L. Note that in each plot the 
last bin is an overflow bin. Here "Other EW" refers to all SM 
backgrounds other than Z/"/* tt and background arising 
from jet^ r misidentification. 



mass" myis , derived from the sum of the observed lepton 
four-momenta and the missing transverse energy. The 
observed distribution of this quantity is dominated by 
the effects of the missing neutrino energy in the tau de- 
cays and experimental resolution. Figure [T] shows the 
observed and fit distributions for the search channels, in- 
cluding the contribution from a Higgs boson signal as 
described below. 

Various uncertainties limit the sensitivity of our search. 
The one with the largest effect is due to the imprecisely 
known tau energy. The distribution of the observed 
transverse momentum of the t m W tv events con- 
strains the ratio of the reconstructed tau energy in the 
observed events to that in the simulation to less than 
1%, but the residual uncertainty allows for shifts in the 
background rrims spectra mimicking a Higgs boson signal, 
particularly for the lower masses considered (m^ < 140 
GeV/c^). At larger Higgs boson masses the search sen- 
sitivity is limited more by other systematic effects con- 
sidered, including the lepton trigger and identification 
uncertainties (2.4% for electrons, 2.7% for muons, and 
4.2% for hadronically decaying taus), the uncertainty in 
the integrated luminosity (6%), Z t^t^ cross sec- 
tion (2.2%), and Higgs boson production cross section 

(5.7%) Hi. 

We represent all the systematic uncertainties by 
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lihood, and eliminate these parameters by maximizing 
the hkeUhood with respect to them. This procedure 
is numerically nearly identical to eliminating them by 
Bayesian marginalization with a Gaussian prior density, 
which takes much longer to compute. 

The resulting likelihood is calculated as a function of 
the Higgs boson signal cross section times branching ratio 
to tau pairs a ■ B, and then converted to a posterior 
density in cr ■ _B assuming a uniform prior density. We 
exclude with 95% C.L. any <j ■ B above which 5% of the 
posterior probability lies. 

The likelihood as a function oi a-B reveals no evidence 
for the presence of a Higgs boson signal, and all nuisance 
parameters remain consistent with their nominal values. 
Figure [T] and other kinematic distributions not shown in 
this Letter all reveal excellent agreement of the observed 
distributions with the predictions. We therefore proceed 
to use the null result to infer upper limits on the Higgs 
boson production cross section. 

Table HIl lists the observed 95% C.L. upper limits on 
a ■ B, and the median upper limits expected under the 
null hypothesis. Figure 2 depicts the results, including 
±1(7 and ±2(7 ranges for the expected limits. We note 
that these bounds would apply to any scalar with similar 
production kinematics decaying to tau pairs. The Higgs 
boson signal distribution shown in Fig. [T] corresponds to 
that excluded at 95% C.L. for tua = 140 GeV/c^. 

We can interpret the upper limits on cr-i? in the context 
of the MSSM parameters tan/3 and uia- The resulting 
excluded regions are shown in Fig. [3] for various assump- 
tions about the sign of the higgsino mass parameteru and 
two extremes for the nature of scalar top mixing [22| , de- 
noted m™°^ and "no mixing." The excluded regions are 
the most stringent published to date in the high tan (3 
region, and are remarkably insensitive to changes in the- 
oretical assumptions due to cancellation of effects in the 
Higgs boson production and decay f23j . 

In summary, we have used a sample of data from the 
Tevatron collider recorded by the CDF II detector corre- 
sponding to 1.8 fb^^ of integrated luminosity to search 
for Higgs bosons predicted in two-Higgs-doublet models, 
via the Higgs boson decays to tau lepton pairs. No ev- 
idence for a Higgs boson signal is observed, and we use 
the null result to infer cross sections excluded at the 95% 
C.L. as a function of the Higgs mass, and 95% C.L. ex- 
cluded regions of the MSSM parameter space tan [3 versus 
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FIG. 2: Observed 95% C.L. upper limits on the cross sec- 
tion for — h/A/H production as a function of niA- The 
grey bands show the median expected limit under the null 
hypothesis, and indicate the ±1- and ±2-standard-deviation 
ranges. 
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FIG. 3: Regions in the MSSM plane of tan /? versus niA ex- 
cluded at 95% C.L. assuming heavy (~1 TeV/c^) sfermions. 
The top panel shows excluded regions for higgsino mass pa- 
rameter /i > 0, and the bottom panel shows excluded regions 
for ^ < 0. Each panel shows the slightly different excluded 
regions for two scalar top mixing scenarios. The solid and 
dashed curves show the previously published bounds, [g] 
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